Decay of Correlations in Fermi Systems at Non-zero Temperature 
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The locality of correlation functions is considered for Fermi systems at non-zero temperature. We 
show that for all short-range, lattice Hamiltonians, the correlation function of any two fermionic 
operators decays exponentially with a correlation length which is of order the inverse temperature 
for small temperature. We discuss applications to numerical simulation of quantum systems at 
non-zero temperature. 
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The locality of correlation functions is a fundamental 
property of a quantum system, with great practical im- 
portance in density functional theory simulation of quan- 
tum systems [1]: if the correlation functions are short- 
ranged, then one can develop fast methods 0{N) for sim- 
ulating a system by solving it in local regions [2,3]. How, 
though, can we know a priori that the correlation func- 
tions will be short-ranged? For quantum systems, it was 
proven recently [4] that the exponential decay in space 
of connected correlation functions follows from the exis- 
tence of a gap between the ground and first excited states. 
If instead the gap vanishes, as occurs at a quantum phase 
transition [5], then long-range correlations may occur at 
zero temperature, while slow fluctuations of the order 
parameters may give rise to exotic physical properties. 

However, another important possibility in which cer- 
tain correlations will decay exponentially in space is to 
consider a system at non-zero temperature. In this paper, 
we prove that correlation functions of fermionic opera- 
tors are short-ranged at non-zero temperature, and we 
show that the correlation length is bounded, for small 
temperature, by a quantity which scales inversely with 
the temperature. The results are general, and apply to 
any lattice Hamiltonian Ti. with finite range R as defined 
below. There are no periodicity requirements and no re- 
quirements of free or weakly interacting particles. 

Previously, it had been suggested that for non- 
interacting particles the correlation length would be in- 
versely proportional to the square-root of the tempera- 
ture [6]. However, later it was argued that, for certain 
specific systems with free particles in periodic potentials, 
the correlation length may be much marger, scaling in- 
versely with the temperature [7]. Thus, the results in [7] 
provide an example showing in this paper we have, in 
fact, derived the best possible bound on the scaling of 
the correlation length with temperature for small tem- 
perature. 

It is important that we consider fermionic operators; 
in contrast, bosonic operators may develop long-range 
order at non-zero temperature. To define what is meant 
by fermionic and bosonic operators, we first define, for 
each lattice site, various fermionic and bosonic opera- 
tors which act only on that site. Possible examples of 



fermionic operators on a site i include electron creation 
and annihilation operators, ipj^^i; examples of bosonic 
operators may include spin operators, phonon operators, 
etc... The defining property of an operator on a single 
site being fermionic or bosonic is that operators on dif- 
ferent sites anticommute with each other if they are both 
fermionic, while they commute with each other if at least 
one of the operators is bosonic. In general, we define an 
operator A to be fermionic if A can be written as a sum 
of products of single-site operators, where each product 
includes an odd number of fermionic single-site operators. 
Thus, ipi and ipiip]ipk are all allowed fermionic operators. 
We define the correlation function at inverse tem- 
perature P by {AB)fj = Z-^Tt[AB exp{-f3H)], where 
Z = Tr[exp(— /3iJ)]. The trace is taken in the grand 
canonical ensemble. Then, we show that if A and B are 
fermionic operators separated by distance I then 



(AS)/3<c||^||||i?||exp(-//0, 



(1) 



where c is a constant and, for large /3, the correlation 
length ^ is of order v(3, with v a characteristic velocity 
of the system. Here, we define the distance I between 
operators to be the shortest distance between any pair of 
sites i,j where some operator on site i appears in A and 
some operator on j appears in B. 

The finite range condition is defined following [4,9]: 
we require that we can write Ti. = ^^ Tii, where i ranges 
over lattice sites, and where the commutator [Tii, O] — 
for any operator O which only acts on sites j which are 
more than a distance R from site i, and where, for some 
constant J, ||7ii|| < J for all i. These conditions are suffi- 
cient to enforce the requirement of a finite group velocity 
below. As an example of the finite range, a hopping term 
4'iipj in the Hamiltonian Tii is allowed if site j is within 
range R of site i . 

To prove Eq. (1), we first express the correlation func- 
tion in terms of anti-commutators using an integral rep- 
resentation. Then, we bound the anti-commutators using 
the finite group velocity [8]. 

Integral Representation of Correlation Function — Let 
A have matrix elements Aij in a basis of eigenvectors of 
Ti, where states i, j have energies Ei, Ej respectively. Let 
{A^)ij = AijS{Ei - Ej - to). Note that A = J dtoAu,. 



The expectation value {A^B)fs — Z ^J2ij^i-^i " 
Ej - u})AijBj,exp{-l3Ei). Similarly, (sA^)^ = 
Z'^ E^ J HE, - E.J - uj)Bj,A,j eicp{~l3Ej). However, 
5{Ei - Ej - uj) exp{-PEj) = 5{E, - Ej - lu) cxp{-/3E,) 
exp(/3cj). Thus, {BA^)p = (A^B}^cxp(/3tj). Therefore, 



{A^B)(: 



1 + exp(/3w 



■{{A^,B})p. 



(2) 



It was essential to use the grand canonical ensemble to 
derive this relation: if instead we worked in the canoni- 
cal ensemble with a fixed number, N, of particles, then 
in the definition of {Ai^B)f3 as a sum over states i,j we 
would require that state i had N particles while in the 
definition of (BAt^)/} we would instead require that state 
j had N particles. If the operator A changes the num- 
ber of particles we would not be able to relate these two 
expectation values as simply as was done above. 

Next^we use (1 +exp(/3w))-i = 1/2-/3-1 EnoddC'^ " 
m7r//3)-i, where the sum ranges over all positive and 
negative odd n. For n > 0, we have (lu — imr/P)~^ = 
i L dtexj>[—{iLu + m:/(3)t\. Similarly, for n < 0, we have 
{uj — imr/ (3)^^ — — «/q dtexp[(ia; + nir/ f3)t\. Thus, 



i, t' > 0; and (3) there exists a constant ci such that 
g{cil,l) is exponentially decaying in / for large I with 
some correlation length ^c- Defining v — Ci , we can 
view w as a characteristic velocity of the system. 

While in [9,4] only bosonic operators were considered, 
for fermionic operators, a similar bound can be proven: 
\\{A{t),B}\\ < \\A\\\\B\\J2j9it,l3)- This proof can be 
proven following the same steps as in [4] with commuta- 
tors replaced by anti-commutators throughout. 

We now combine the bound on anti-commutators with 
Eq. (5). The physics is as follows: Eq. (4) gives the 
correlation function as an integral over time. However, 
for times t large compared to (3, the integral is cut off 
by the exponential; however, for times t of order /3 or 
less, the anti-commutator of Eq. (4) is small for I large 
compared to vp. 

More precisely, note that since A, B are separated by 
a distance /, then {A, B} = 0. Thus, 



cil 



'Pjo 
, i 



dt 



1 



1 -I- exp(/3w) 



exp(— 7rt//3) 
1 - exp(-27ri//3) 
cxp(— 7ri//3) 



dt 



1/2+ (3) f^Jcii 1 - cxp(-27ri//3) 



\{AB)p\< (6) 
{{A{t)-A{-t),B})^\ + 

{{Ait)-A{-t),B})^\. 



/•OO 

y / dt exp(— n7rt//3)(cxp(ia;i) — cxp(— iwt)), 



nodd,n>0 



where now the sum ranges only over positive, odd n. 

Now, we define the time evolution of operators 
by A{t) = cxp{im)Aexp{-int). Thus, A^{t) = 
exj>{iLut)A^. Combining Eqs. (2,3), we get 

{A^B)^^1/2{{A^,B})(,+ (4) 

-B rdt /f--;/^) ({A.(t)-A.M),B}),, 
P Jo 1 - cxp(-27rt//3) 

where we have used I]„odd,ri>o'^^P(^"^V/^) ^ 
exp(— 7rt//3)/(l — exp(— 27rt//3)). Finally, we integrate 
Eq. (4) over iv to get 

{AB)f,^l/2{{A,B})p+ (5) 

P Jo 1 - exp(-27ri/p) 

Finite Group Velocity — In [9,4], it was proven that if 
we have a Hamiltonian obeying the finite range condition 
above, then for any two bosonic operators A, B, separated 
in space by a distance I, one can bound the operator norm 
of the commutator: ||[^(i),B]|| < \\A\\\\B\\J2j g{t,lj), 
where the sum ranges over sites j which appear in op- 
erator B and Ij is the distance from j to the closest 
site i in the operator A. Here g{t,lj) is a function with 
the following properties: (1) g is symmetric in t so that 
g{t,l) = gi-t,l); (2) g^tj) < {t/t')g{t' ,1) for t < t' and 



Note that \{{A{t) - A{-t),B})f3\ < 4||A||||B||. Thus, 
the integral over times t > ciZ in Eq. (6) is bounded 
in absolute value by 4||A||||B||/3-i /^'^ exp(-7rt//3)/(l - 

exp{-27Tt/(3)) < 4||yl||||B||^-icxp(-^ciV/3)/(l - 
exp(-27rci//;9)). 

For 
times t < cil, we use the bound on the anti-commutator: 
\{{A{t),B})p\ < \\{A{t),B}\\ < [t/icj)]j:,gicil,l,)- 
Also, exp(-7rt//3)/(l - exp(-27rt//3)) < /3/(27ri). Thus, 
the integral over times t < cil in Eq. (6) is bounded in ab- 
solute value by \\A\\\\B\\j:^ g{cj,l,) j;;'' dt{2TTCil)-^ = 
{2n)-^\\A\\\\B\\j:^g{c,l,l,). 

Therefore, 



\{AB)^\<\\A\\\\B\ 



Ej gjcii, ij 

2tt 



4 exp(— 7rci^//3) 
TT 1 — exp(— 27rci///3) 



Since Ij > I, g{cil,lj) decays exponentially for large 
I as exp(— Z/^c")- Also, exp(— ttciZ//?) decays exponen- 
tially with correlation length u/S/tt. Thus, the correlation 
length ^ is bounded by the maximum of ^c and w/J/tt, 
so that for small temperature ^ < vP/tt. Thus, Eq. (1) 
follows. 

It is interesting to compare Eq. (7) to the result found 
in [4] for systems in the ground state with a gap /S.E 
between the ground and first excited states. There, the 
resulting correlation length was the maximum of ^c" ^^nd 
2v/l^E. 

Importance of Lattice Structure — This bound was de- 
rived for systems with a lattice structure. To illustrate 



(7) 



the importance of a lattice, consider free fermions on a 
lattice in one dimension, with Hamiltonian til^jipi^i + h.c. 
The Fermi velocity depends on filling, and is maxi- 
mum at half-filling. Thus, the Fermi velocity can be 
bounded, and is at most of order t. Now, consider free 
fermions in free space one dimension, with Hamiltonian 
ip'^ {x)[d'^ / {2m)]il;{x) . The Fermi velocity again depends 
on the particle density, but can be increased without limit 
as the particle density is increased. 

This is the physical reason why, for systems which 
are not on a lattice, one cannot (without knowing more 
about the density and other details of the system) bound 
the anti-commutator. One cannot, from the Hamiltonian 
alone, provide a velocity v such that the anti-commutator 
{A{t), B} is small for t < I /v. 

Discussion — We have proven a bound on the corre- 
lation of fermionic operators on non-zero temperature, 
Eq. (7) . There are specific examples [7] which show that 
this bound on the scaling of the correlation length with 
temperature is the best possible (although the prefactor 
w/tt in the correlation length ^ — vP/tt might not be the 
best possible). Similarly, while it had been shown that 
for periodic, non-interacting insulators in one-dimension 
the correlation length scales inversely with the square- 
root of the gap [10], there are again specific examples of 
non-interacting systems [7] which show that the bound 
[4] on the scaling of the correlation length with the gap 
is also the best possible. The case of the canonical, as 
opposed to grand canonical, ensemble is an interesting 
future problem. 

In addition to the basic interest in this result in quan- 
tum statistical mechanics, this result is of importance 
in quantum simulation using density functional theory. 
Also, the integral representation, Eq. (5) may prove use- 
ful as a means of computing the density matrix for free 
fermions at non-zero temperature; for these systems, if 
A{Q) is the fermionic creation operator on a single site 
and _B is a fermionic annihilation operator, then A{t) 
is a sum of fermionic creation operators and the anti- 
commutator is trivial to compute. This is a matter for 
future research. 
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